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Abstract - The free piston energy converter (FPEC) is a new 
type of power unit for hybrid electric vehicle. In a FPEC, an 
internal combustion engine driven permanent magnet linear 
generator (PMLG) may be utilized to produce electrical energy. 
This paper discusses the issues that are pertinent to the design 
of a linear permanent magnet generator for application in a 
FPEC. To achieve the required high power density and low 
moving mass, a tubular machine equipped with modular 
windings and axially magnetized is considered to be the most 
appropriate. Finite Element Analysis are adopted to analyse the 
PMLG. Magnetic field distribution, EMF, thrust force and No- 
load corl loss are obtained through the finite element method. 
An effective tool for design optimization is provided. It is shown 
that the machine design can be optimized with respect to three 
key dimensional ratios while satisfying other performance 
specifications. 


I. INTRODUCTION 


Efficiency improvements in order to reduce the production 
of CO2are required in all technical areas of life to limit the 
anticipated climate changes to an extent sustainable for 
humanity [1]. Current approaches in the field of 
transportation cover both increasing the efficiency of 
conventional drives and the development of hybrid, battery, 
fuel cell and range extender concepts, and combinations of 
these [2]-[5]. This requirement is met particularly well by 
Free piston energy converter (FPEC) [2],[5],[6]. It is capable 
of transforming chemical energy into electrical energy by 
means of a combustion process. This is achieved by keeping 
the system frequency constant and adapting to the power 
demanded by variation of the stroke and compression. 
Furthermore, these system characteristics mean that the free- 
piston linear generator can be operated with both 
conventional fuels such as petrol, diesel and gas and with 
alternative fuels such as sun fuel, synthetic fuel, hydrogen 
etc[7]-[9]. 
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Fig. 1 Schematic diagram of the FPEC 


The schematic diagram of the FPEC is shown in Fig. 1. 
The two ends of the translator of a linear generator are 
connected to the two pistons of two oppositely placed 
combustion chambers, respectively. The reciprocating 


ignition and compression processes in the two chambers 
drive the translator of the linear generator to move backward 
and forward, thus, electric energy is produced in the 
generator, which will be supplied to a motor, after 
rectification and inverting, to drive the vehicle. As will be 
evident from Fig. 1, Thus,The linear generator is an 
important part in the FPEC, and in this paper a tubular 
axially magnetized magnet linear generator scheme is 
investigated. 


Il. STRUCTURE OF THE TUBULAR AXIAL FLUX PM LINEAR 
GENERATOR FOR FPEC 


A 3-phase, 12-slot/11-pole tubular linear generator 
equipped with a modular stator winding in which the coils of 
each phase are disposed adjacent to each other, as shown in 
Fig. 2, has been identified as being the most appropriate for 
this application. It has two parts: stator and translator. The 
stator is composed of the stator lamination and winding, and 
the translator consists of axially magnetized permanent 
magnets, iron segments, and stainless steel. 

silicon steel is the best material for low reluctance 
magnetic circuit material, so the stator is constructed from 
non oriented silicon steel laminations of grades 50W470. 
The typical specifications of 50W470 is shown in Table 1. 
The permanent magnets are made of a high flux density rare 
earth Neodymium Iron Boron with a remanence of 1.21T and 
a relative recoil permeability of 1.06. the iron segments is 
mild steel 45 which are sandwiched between the permanent 
magnets to get better magnetic performance. The stainless 
steel material which has a relative permeability around 1 is 
chosen as the shaft material. The low relative penneability 
shaft causes almost all fluxes to flow through the stator 
core.To decrease the total weight and increase the specific 
power, the long translator structure is adopted in this paper. 


Phase B Phase C 


Phase A 
Permanent 
magnet 


winding 
Stator lamination 


Fig. 2. Schematic of 3-phase, 12-slot/11-pole, axially magnetized magnet 
linear generator. 
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Table 1 
Lamination material property 
i Max.core loss Max. Inductio 
Grade hk pa (W/kg) nat 
m) gcm) LOT LST 5,000A/m 
50W470 0.5 Ta 2.00 4.70 1.65 


A. Armature reaction 

The magneto-motive force (mmf) distribution in the 
slotted stator of a tubular PM machine may be represented 
using a current model, assuming that the yoke and teeth are 
infinitely permeable. According to Ampere’s law, Figs. 3(a) 
and Figs. 4(a) show the actual MMF distribution of the three 
phase winding of the 3-phase, 12-slot, 11-pole modular 
machine with single-layer and double-layer windings, 
respectively. Harmonic contents of the three phase winding 
of 12-slot, 11-pole modular machine with single-layer and 
double-layer windings are provided in the Figs. 3(b) and 
Figs. 4(b). Fig. 3(b) shows a worsening of the harmonic 
contents when the winding is transformed into a single-layer 
one. In particular, the increase of the subharmonic of order v 
= | is evident. The lower and higher order mmf harmonics 
travel at different speeds to the armature, and induce eddy 
current losses in both the PMs, iron, and stainless steel, and 
may cause other undesirable effects such as localized core 
saturation, noise, and vibration. 
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(a) (b) 
Fig. 3. MMF distribution of a 12-slot 11-pole linear generator with single 
layer. (a) Actual waveform. (b) harmonic content 
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Fig. 4. MMF distribution of a 12-slot 11-pole linear generator with double 
layer. (a) Actual waveform. (b) harmonic content 


B. Flux Distribution 

Calculate the induced voltage in the coils and map the 
flux density in the generator parts by flux calculation. The 
voltage induced in the coils can be derived from the flux 
changing with respect to time. The shape and the dimension 
of teeth, the back iron of the stator, and air gap are 


determined by the flux density. The magnetic potential 


vector A , instead of flux density, is often used in the field 
analysis. The differential operation to this vector potential 
gives: 


> > 
V-B=V-(Vx 4) =0 (1) 
For problems considering saturable material with 


permanent magnets, the constitutive relation for the magnetic 
fields is: 


as 


B=yH+ uM, (2) 
MyM, = B, (3) 


where [Up a is absolute permeability of free space, M, is 


remanent intrinsic magnetization vector, and B, is remanent 
flux density. The flux flows through an area can be 
calculated as: 


y= [B-nds (4) 


C. Output Voltages 
Electro-motive force (EMF) of one phase is calculated 
from the linkage flux in the coils by Faraday's Law. 


d 
e,=-NY (5) 
dt 
where N is the number of turns. 
The terminal voltage U, for a phase can be expressed as: 


Us =e, — Rsis + L, di, (6) 
` dt 
where R, is phase resistance, i, is the phase current,, L, is 
phase inductance. 
Fig. 5 shows the equivalent per phase circuit diagram of 
the PMLG. 


Fig. 5. Equivalent per phase circuit diagram 


D. Output Power and Efficiency 
Input power is defined as follows if the mechanical loss is 
also neglected, 


Fn = Eo + Pret Pou + Pead (7) 
Where P „is the output power. p,p, is the core loss, which 


can be calculated by using empirical formulas , p,, is copper 
loss which is proportional to the square of load current and 
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Paa iS eddy-current loss. 
The output power is calculated as following, 


P „=U, +U, M, +UA, (8) 


Where P „is the average output power, U,, U, and U, is 
voltage RMS of phase A, B and C, respectively, Z, J, and Z, 
are current RMS of phase A, B and C respectively. 

The efficiency of a PMLG is obtained as follows, 


y = 2x100% (9) 


in 


E. Cogging Force 

The cogging force is due to the magnetic attraction 
between the permanent magnets mounted onto the shaft and 
the stator teeth. The force attempts to maintain the alignment 
between the permanent magnet and the teeth. It becomes an 
important parameter since its peak value is significantly 
high. Virtual work method is presented for computing the 
cogging force in the x direction, expressed as follows, 


F =-— (10) 


Where W is the surface energy of mover on air-gap, 


expressed as, 
k = > 
w= || f (#-a8) fs, (11) 


i=l g§ 
Where k is the total number of local triangular elements on 
mover surface, H , is the magnetic field intensity, and B, is 
the magnetic flux density at element i. 


Ul. DESIGN OPTIMIZATION 


The design of the machine can be optimized with respect 
to the design parameters shown in Fig. 6. For a given stator 
outer radius Ro, The generator design parameters, such as the 
mover inner radius R;, the dimensional ratios T,/T,, the 
magnet thickness Am, turns of coils N., and the airgap length 
g are determined iteratively considering the overall generator 
size loss, efficiency and other parameters to arrive at an 
optimal system design. 
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Fig. 6. Design parameters of modular tubular linear generator 


In general, the air-gap length g is assumed to be constant 
at 1 mm, since although a smaller air-gap length would also 
improve the performance, it is limited by manufacturing 
tolerances. the performance improves as h,, is increased. 
However, an increase in the volume of rare-earth PM 
material increases the cost and results in a heavier armature, 
which is usually undesirable for amoving-magnet machine. 
In this study, therefore, the magnet thickness is fixed at 5 
mm to produce an acceptable air-gap flux density and thrust 
force. Fig. 7 shows the variation of thrust force as a function 
of Imax and N, assuming R=0.029m. Fig. 8 shows the 
variation of thrust force as a function of /,,,, and R;. As can 
be seen, irrespective of the Fig. 7 or Fig. 8, within specified 
limits, as three phase current, the number of turns or the 
mover inner radius increases, the thrust force almost 
increases with linearly. 
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Fig. 7. thrust force as a function of Jmax and Ne. 
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Fig. 8. thrust force as a function of Jmax and Ri. 


Fig. 9 and Fig. 10 show phase-to-neutral back-EMF 
waveforms and the rms value of induced voltages for the 12- 
slot 11-pole linear generator, respectively. it is possible to 
observe that this EMF exists an optimal ratio of 1,/t,°0.7 
which yields the maximum rms value of induced voltages 
and the best sine wave. The no load core loss versus the 
variation of the dimensional ratios T,,/T, is shown in Fig. 11. 
As can be seen, the no load core loss is reduced gradually 
with the Tm/Tp increased. 
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Fig. 9. Phase-to-neutral back-EMF waveforms for the 12-slot 11-pole linear 


generator calculated by FEA 
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same conclusion can be obtained, and only the current phase 


changes. 
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Fig. 12. Output power and efficiency for PMLG at 10 m/s 
IV. FINITE ELEMENT ANALYSIS 


A linear generator design are presented here to illustrate 
the application of the finite element method. Fig. 13(a) and 
(b) show, respectively, the no-load and full-load flux 
distributions of the machine, in which end-effects associated 
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Fig. 10. Variation of back-EMFs with dimensional ratios Tm/ Tp 
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Fig. 11. Variation of no-load core loss with dimensional ratios Til Tp 


The output power and the efficiency for PMLG in 
different loads are shown in Fig. 12. As can be seen that the 
output power and efficiency climb up and then decreases 
with an increasing load resistance. When the load resistance 
is 2 Q, 8.7kW output power is got and the efficiency reached 
93.96% at a constant speed of 10 m/s. For RL loads, the 


with the finite stator length are considered. 


B[tesla] 


1. 8763e+000 
1. 7423e+000 
1. 6083e+000 
1.4742e+000 
1, 3402e+000 
1, 2062e+000 
1.0722e+000 
9, 38156-0081 
8, 4132-001 
6. 7011-001 
5, 3609e-001 
4, 82072-0014 
2. 6805-001 
1. 3403-001 
5, 99306-0065 


(a) No load 
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(b) full load 
Fig. 13. Flux distributions in 3-phase, a 12-slot 1 1-pole linear generator. 
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Fig. 14 shows the Br and Bz components for the no load 
magnetic flux density of the linear generator. Induced 
electromotive force when load is 2 and the speed is 10m/s is 
shown in Fig. 15. It is clear that the three phase voltage is 
quite symmetrical. 
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Fig. 14. B, and B, components for the no load magnetic flux density of the 
linear generator 
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Fig. 16. Thrust force waveforms with sinusoidal current excitation at 10m/s. 


Fig. 16 shows the thrust force waveforms resulting from 
sinusoidal current excitation and the average thrust force is 
2.55 KN at 7 „=80A. The peak-to-peak thrust force is in the 
vicinity of 0.21 kN at J,,=0A, corresponding to 
approximately 8.2% of the linear rated thrust force. A variety 
of techniques are available to further reduce the thrust force 
ripple. For example, the magnet span can be adjusted in 
order to nearly eliminate the lowest frequency component of 
the thrust force, Bifurcated slots dummy slots in stator 


teeth. 


V, CONCLUSION 


The suitability of employing a tubular PM linear generator 
equipped with axially magnetized magnets in a free piston 
energy converter has been discussed. It has been shown that 
this machine topology has attractive features in terms of its 
high power density, high efficiency and low cogging force, 
which are essential for achieving a high output power from 
an FPEC. An effective tool for design optimization has been 
provided. Magnetic field distribution, EMF, thrust force and 
No-load core loss has been obtained through the finite 
element method. The output power of the PMLG has been 
calculated in different loads, and the efficiency reached 


93.96% at 2 Q condition. It is shown that the machine design 


can be optimized with respect to three key dimensional ratios 
while satisfying other performance specifications. 
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